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for  the  Actual  Cases. 

1  •   P ■■  ^ ^ pcse  of  the  present  report  and  suramary  of  its  re s ult s  . 

In  Reoort  I  a  mathematical  basis  vras  laid  for  the 
numerical  solution  of  flow  problems  in  rivers.   The  ultimate 
aim  is  to  carry  out  numerically  the  solution  of  a  flood  problem 
for  the  CtLio  River  in  a  concrete  case  'oj   using  an  appropriate 
digital  comDuter,  to  compare  the  results  vrith  the  observations, 
and,  in  general,  to  study  the  feasibility  of  such  methods  of 
attacking  this  type  of  problem;  in  addition,  the  problem  of 
floods  at  the  Junction  of  the  Ohio  and  the  Mississippi,  and 
problems  coneerninf  the  regulation  of  the  Tennessee  River 
through  controls  at  the  Kentucky  Dam  are  to  be  solved 
n\americally.   In  all  of  these  cases  it  is  necessary  to  make 
use  of  a  considerable  bulk  of  observational  data--cross- 
sections  and  slones  of  the  channels,  measurements  of  river 
depths  and  discharges  as  functions  of  time  and  distance  down 
the  river,  drainage  areas,  observed  flows  from  tributaries, 
etc. --in  order  to  obtain  the  information  necessary  to  fix  the 
coefficients  of  the  differential  equations  derived  in  Rei-^ort  I 
and  to  fix  the  initial  and  boundary  conditions.   This  in  itself 
is  a  task  with  some  complexities,  and  it  also  tal:es  time  and 
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the  coor)eration  of  several  oroups  of  people.   In  the  nieantiTiG, 
therefore,  it  was  thoug-ht  wise  to  try  out  the  numerical  inethods 
proposed  in  Report  I  on  three  problems  which  are  simplified 
versions  of  actual  problems,  and  to  carry  out  the  solutions 
by  usin;-^  ordinary  hand  calculators.   The  present  Rex^ort  II 
has  as  its  purpose  the  presentation  of  the  solutions  in  these 
special  cases,  together  with  an  analysis  of  their  bearing  on 
the  concrete  problems  for  actual  rivers. 

The  models  chosen  correspond  in  a  rough  general  way 
(a)  to  two  tyces  of  flow  for  the  Ohio  and  (b)  to  the  Ohio  and 
^'lississiopi  at  their  junction.   Rivers  of  constant  slope,  with 
rectangular  cross-sections  having  a  uniform  breadth  and 
roughness  coefficient   are  assumed.   In  this  way  differential 
equations  with  constant  coefficients  result.   The  values  of 
these  quantities  are,  however,  taken  to  correspond  in  order  of 
magnitude  with  those  for  the  actual  rivers.   In  the  model  of 
the  Ohio,  for  example,  the  slope  of  the  channel  was  assumed 
to  be  0,^   ft. /mile,  the  quantity  n  (the  roughness  coefficient  in 
Manning's   formula)  was  given  the  value  0.03,  and  the  breadth 
of  the  river  was  taken  as  1000  feet.   (These  values  were 
recommended  to  us  by  the  engineers  on  the  basis  of  their 
knowledge  of  the  Ohio.)   It  was  assumed  that  a  steady  lUiiform 
flow  with  a  depth  of  20  ft.  existed  at  the  initial  instant 
t  =  0,  and  that  for  t  >  0  the  depth  of  the  water  x-j&s    increased 
at  the  point  x  =  0  from  20  ft.  to  ItO  ft.  viithin  L|.  hours  and 
then  held  fixed  at  the  latter  value.   The  problem  was  to 
determine  the  flow  doi«jnstream,  i.e.  the  depth  y  and  flow 
velocity  v  as  functions  of  x  (for  x>  0)  and  t.   The  resulting 
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stage  and  discharge  profiles  are  given  for  several  times  in 
Figures  1  and  2.   In  a  similar  fashion  the  model  for  the 
junction  of  the  Ohio  and  Mississippi  was  set  up.   (Pigxires  l^ 
and  5  display  the  stage  and  discharge  profiles  for  various 
times . ) 

Much  valuable  Information  and  insight  was  gained  from 
the  study  of  these  models,  especially  from  the  model  of  the 
Ohio.   This  will  be  discussed  in  detail  later  on,  but  a  number 
of  observations  mi[^ht  be  made  at  this  point.   In  tlie  first 
place,  the  fact  that  the  solution  of  the  rlraplified  flood 
problem  in  the  Ohio  could  be  carried  out  numerically  by  hand 
computation  over  a  considerable  range  of  values  of  the 
distance  and  time  (values  at  900  net  points  in  the  x,t-plane 
were  determined  by   finite  differences)  shotirs  that  the  problems 
are  well  within  the  capacity  of  modern  calculating  equipment. 
It  mi<^ht  be  added  that  the  special  case  chosen  for  a  flood 
in  the  Ohio  was  one  in  which  the  rate  of  rise  at  the  starting 
point  upstream  was  extremely  high  (5  feet  per  hour,  in 
comparison  with  the  rate  of  i-ise  dxirlng  the  big  flood  of 
19^4-5  Hhich  was  never  larger  than  0.7  feet  per  hoiar  at 
V/heeling,  West  Virginia),  so  tiiat  a  rather  severe  test  of 
the  finite  difference  method  was  made,  in  view  of  the  rapid 
changes  of  the  basic  quantities  in  space  and  time. 

The  decisive  point  in  estimating  the  magnitude  of  the 
computational  work  In  using  finite  differences  is  the  number 
of  net  points  needed,  and  our  model  indicates  that  an  interval 
/ \ X  not  smaller  than  10  miles  along  the  river  and  an  interval 
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/\ t  not  less  than  C.3  hours  in  time  in  a  rectangular  net  in 
the  x,t-plane  will  yield  results  that  are  sufficiently 
accurate.   (Of  course,  the  actual  problems  for  the  Ohio  vjill 
involve  empirical  coefficients  in  the  differential  equations 
and  other  empirical  data,  which  will  have  to  be  coded  for 
calculating  machines,  but  this  will  have  no  ,e;reat  effect  on 
these  estimates  for  /\  x  and  iniiP-ht  under  extreme  flood  condi- 
tions reduce  /\t     by  a  factor  of  .$.)      In  our  model  calcula- 
tions we  used  a  uniform  space  interval  /\x  ,  alonp  the  river's 
length,  but  we  are  investigating  the  feasibility  of  using  the 
location  of  the  gaging  stations  as  net  points.   In  our  model 
calculations,  we  exneriment  with  various  finite  difference 
methods,  interval  sizes  and  analytic  approximations  in  order 
to  get  the  practical  experience  necessary  to  select  the 
technique  to  be  used  for  automatic  iiachine  computation,   '^'e 
plan  to  use  our  ideas  in  some  trial  desk  calculations  made  vjith 
the  physical  data  fo^^  the  Ohio  and  the  19u-5  flood  (the  basic 
data  have  been  provided  by  the  Corps  of  Engineers  and  is  now 
being  reduced  at  ".Tew  York  University  to  a  form  suitable  for 
computation) . 

A  second  computation'  for  the  Ohio  mr^del  was  performed 
for  a  steady  progressing  wave  {hO   ft.  stage  upstream  and  20  ft. 
stage  downstream) ,   This  work  indicated  that  an  interval  ,A^x 
of  10  miles  yielded  a  stage  prediction  accurate  to  within  0.6 
Tjer  cent  in  a  7  hour  forecast. 


To  provide  us  with  a  check  on  the  accuracy  of  our 
computational  methods,  in  a  problem  vmich  could  be  solved 
analytically  (see  sec.  2.}.\.   and  App.  II). 
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As  a  model  fov  the  problem  of  calculating:  what  happens 
at  the  junction  of  two  major  streams,  vre  worked  with  three 
separate  river  stretches  (e.g.  the  Ohio,  the  upstream  side  of 
the  Mississippi,'  the  downstream  side  of  the  Mississippi).   For 
each  region  the  same  finite  difference  scheme  was  used  as  in 
the  model  of  the  Ohio.   In  addition,  it  was  necessary  to 
determine  the  values  of  stage  and  discharge  at  the  common 
junction  point  from  the  knowledge  that  the  levels  of  the  three 
branches  agree  there  and  that  the  water  which  flows  into  the 
junction  also  flows  out. 

A   much  better  understanding  was  gained  from  the  calcula- 
tions for  the  model  o.';'  the  Ohio   of  the  relation  between  the 
methods  used  by  the  engineers  in  the  Ohio  Piver  Division  in 
Cincinnati  (and  other  engineers  as  v/ell)  for  predicting  flood 
stages,  and  the  methods  explained  in  Report  I,  which  rcake  use 
of  the  basic  differential  equations.   At  first  sight  the  two 
methods  seem  to  have  very  little  in  comm.on,  though  both,  in 
the  last  analysis,  must  be  based  on  the  laws  of  conservation 
of  mass  and  momentuin;  indeed,  in  one  important  respect  they 
even  seem  to  be  somewhat  contradictory.   Hie  methods  used  in 
engineering  practice  (which  make  no  direct  use  of  any 
differential  equations)  tacitly  assume  that  a  flood  wave  in  a 
long  river  such  as  the  Ohio  propagates  only  in  the  dovjnstream 
direction,  v.'hile  the  basic  theory  of  the  differential  equations 
we  use  (as  was  explained  in  Report  T)  tells  vis  that  a 
disturbance  at  any  point  in  a  river  flowing  at  sub-critical 
speed  (the  normal  case  in  general  and  al^^Jays  the  case  for  the 


Ohio)  will  propagate  as  a  wave  traveling  upstream  as  well  as 
downstream.   Not  only  that,  the  speed  of  propagation  relative 
to  the  flowing'  stream,  as  defined  by  the  differential  equations, 
is  >/gy  for  small  disturbances  and  this  is  very  much  larger 
than  the  propagation  speed  obtained  by  the  engineers  for  their 
flood  wave  traveling  doi-jnstream.   Tlaere  is,  however,  no  real 
discrepancy.   The  method  used  by  the  engineers  can  be 
Interpreted  as  a  Method  which  yieD.ds  solutions  of  the 
differential  equations,  with  certain  terms  negD.ected,  that  are 
good  approximations  (though  not  under  all  circumstances,  it 
seems)  to  the  actual  solutions  in  some  cases,  among  them  that 
of  flood  x-javes  in  a  river  such  as  the  Ohio.   The  neglect  of 
terms  in  the  differential  equations  in  this  approximate  theory 
is  so  drastic  as  to  make  the  theory  of  characteristics,  from 
which  the  properties  of  the  solutions  of  the  differential 
equations  vrere  derived  in  Report  I.  no  longer  available.   The 
numerical  solution  presented  here  of  the  differential  equations 
for  a  flood  wave  in  a  model  of  the  Ohio  yields,  as  we  have 
said,  a  wave  the  front  of  which  travels  downstream  at  the 
speed  \fgj,    but  the  amplitude  of  this  forerunner  is  quite  small, 
while  the  portion  of  the  wave  with  an  amplitude  in  the  range 
of  practical  Interest  is  found  by  this  method  to  travel  with 
essentially  the  same  speed  as  would  be  determined  by  the 
engineer's  aoi^roximate  method,   ''•'hat  seems  to  happen  is  the 
following:   small  forerunners  of  a  dist\irbance  travel  with  the 
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speed  /gy   relative  to  the  flowing  stream,  but  the  friction 
forces  act  in  such  a  way  as  to  decrease  the  speed  of  the  main 
portion  of  the  disturbance  far  below  the  values  given  by  /"gy, 
i.e.  to  a  value  corresponding  roughly  to  the  speed  of  a  steady 
progressing  wave'   that  travels  unchanged  in  form.   (One  could 
also  interpret  the  engineering  method  as  one   based  on  the 
assumption  that  the  waves  encountered  in  practice  differ  but 
little  from  steady  progressing  waves.) 

This  analysis  of  the  relation  between  the  methods 
proposed  in  ^e^gort    I  and  those  commonly  used  in  engineering 
practice  indicates  why  it  may  be  that  the  latter  methods, 
while  they  furnish  good  results  in  the  most  importaxit  cases, 
fail  to  mirror  the  observed  occurrences  in  other  cases.   Por 
example,  the  problem  of  what  happens  at  a  major  junction,  and 
various  problems  arising  in  connection  with  the  operation  of 
such  a  dam  as  the  Kentucky  Dam  in  the  Tennessee  River  seem  to 
be  cases  in  which  the  engineering  methods  do  not  work  l^rell. 
Tliese  are,  in  all  likelihood,  cases  in  which  the  motions  of 
interest  depart  too  much  from  those  of  steady  progressing  waves, 
and  cases  in  which  the  propagation  of  waves  upstream  is  as  vital 
as  the  propagation  downstream.   Thus  at  a  major  junction  it  is 


As  can  be  seen  from  Report  I,  the  propagation  speed,  Tgy, 
which  is  determined  without  reference  to  the  friction  term, 
would  not  be  affected  by  any  change  in  the  order  of  magnitude 
of  the  friction  force. 

A  steady  progressing  wave  is  a  flow  in  which  the  depth  and 
discharge  are  functions  only  of  (x  -Ut ) ,  with  U  a  constant 
velocity.  Tn.e   propagation  speed  U  is  determined  from  the 
friction  term  (see  sec.  2,li). 


clear  that  considerable  effects  on  the  upstreain  side  of  a  main 
stream  are  to  be  ex-oectod  when  a  large  flow  from  a  tributary 
occurs.   Tn  the  same  way,  a  da:a  in  a  stream  (or  any  obstruction, 
or  change  iii   cross-section,  etc.)  causes  reflection  of  waves 
upstream,  and  ne.s-lect  of  such  reflections  might  well  cause 
serious  errors  on  some  occasions.   On  the  other  hand,  the 
method  of  numerical  solution  based  directly  on  the  differential 
eciuations  is  applicable  in  any  of  these  cases,  and  it  seems 
reasonable,  as  we  propose  doing,  to  test  it  out  as  a  practical 
method.  / 

The  above  general  description  of  what  happens  when  a 
flood  vrave  starts  down  a  long  stream.--in  particular,  that  it 
has  a  lengthy  front  portion  xv'hich  travels  fast,^ut  has  a  small 
amplitude,  while^  the  _main  j)art  of  the  disturbance  moves  much 
more  slox>rly--has  an  important  bearing  on  the  que ct ion  of  the 
proper  approach  to  the  numerical  solution  by  finite  differences. 
It  is,  as  vje  shall  see  shoi'tly,  necessary  to  calculate--or  else 
estimate  in  some  way--the  motion  up  to  the  front  of  the 
disturbance  in  order  to  be  in  a  position  to  calculate  it  at  the 
places  and  times  where  the  disturbances  are  large  enough  to  be 
of  practical  interest,   Tliis  means  that  a  large  number  of  net 
points  in  the  finite  difference  mesh  in  the  x,t-plane  lie  in 
regions  i>rhere  the  solution  is  not  of  great  interest.   Since  the 
fixing  of  the  solution  in  these  regions  costs  as  -nuch  effort  as 
for  the  regions  of  rreater  interest,  the  differential  eqiiation 
method  is  at  a  certain  disadvantage  by  comparison  with  the 
conventional  method  in  such  a  case.   However,  it  is  possible  to 
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determine  analytically  the  character  of  the  front  of  the  wave 
and  thus  estimate  accurately  the  places  and  times  at  which  the 
wave  amplitude  is  so  small  as  to  be  negligible;  these  regions 
can  then  be  regarded  as  belonging  to  the  regions  of  the  x,t- 
plane  where  the  flov/  is  undisturbed,  vrith  a  corresponding 
reduction  in  the  number  of  net  points  at  which  the  solutions 
must  be  calculated.   A  method  which  can  be  used  for  this 
purpose  has  been  devised  by  G.  V-hitham  and  A.  Troesch.  and  a 
description  of  it  is  given  in  Appendix  I  of  this  report.   If  a 
modern  high  speed  digital  computer  were  to  be  used,  however,  it 
would  not  matter  very  much  whether  the  extra  net  points  in  the 
front  portion  of  the  wave  were  to  be  included  or  not:   many 
such  machines  have  ample  capacitj-  to  carry  out  the  necessary 
calculations. 

~n  sec.  2  a  description  of  the  calculations  made  for  the 
model  of  the  Ohio  is  given,  including  a  discussion  of  various 
difficulties  which  occurred  for  the  flood  wave  problem  near 
the  front  of  the  dist-urbance,  and  particularly  at  the  beginning 
of  the  wave  motion  (i.e.  near  x  =  0,  t  =  0).  and  an  enumeration 
of  the  features  of  the  calculation  which  are  certain  to  play  a 
similar  role  in  the  more  complicated  case  of  the  Ohio  as  it 
really  is.   In  sec.  3  we  give  a  description  of  the  method  used 
and  the  calculations  made  for  a  model  of  a  flood  comijig  dOT-.tn 
the  Ohio  and  its  effect  on  passing  into  the  Mississippi.   In 
Appendix  I,  as  was  noted  above,  the  method  of  C-.  ^"hitham  and 
A.  Troesch  for  dealing  with  the  front  oortion  of  a  vjave  traveling 
down  a  river  is  described.   In  Appendix  11,  a  discussion  of 
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steady  flows  and  of  steady  progressing  waves  is  given.   These 
flows  were  used  to  deterraine  convenient  initial  conditions,  to 
test  our  computing  procedu.res  and  to  furnish  the  asymptotic 
behavior  (for  large  time)  of  the  flows. 

Furthermore,  in  order  to  understand  more  full^i^  the 
general  behavior  of  flood  waves  in  open  channels  and  especially 
why  the  speed  fgy  is  not  the  observed  propagation  velocity  for 
normal  flood  problems,  some  theoretical  Investigations  using 
the  concepts  of  steady  and  steady  progressing  flow  were  carried 
out  by  G.  Norikawa,  A.  Troesch,  and  G.  Whitham  and  will  be 
reported  on  later,   'rheir  concliisions  also  have  important 
bearing  on  the  practicability  of  making  flood  forecasts  for 
extended  periods  of  time,  from  the  given  data. 

2,   Model  of  the  Ohio. 

2.1  Formulation  of  the  flood  problem.   We  consider  a  simple 
model  x^fhich  approximates  the  average  characteristics  of  the 
Ohio  River.   The  I'iver  was  assumed  to  have  a  rectangular  cross 
section  1000  ft.  in  width.   The  slope  S  of  the  river  bottom 
was  taken  as  .5  ft,/-'^ile  and  Manning's  form.ula  was  used  for  the 
friction  slope  S^, 
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Hayashi,  T. ,  l^lathematical  theory  of  flood  waves.   Froc.  of  the 
1st.  Japan  Nat'l.  Cong,  for  Appl .  Mech.,1951. 

Deymie',  Fh.,  Fropagation  d'lone  intumescence  allongee  (probleme 
aval).   Int."cong.  for  Appl.  Mech.,5th,  Proceedings, 
Sept.  1938,  N.  Y. ,  ^iiey,  1939. 
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(since  the  hydraulic  radius  R  for  a  rectangular  channel  of 
width  B  and  water  of  depth  y  is 

B 
The  choice 

,  =  ,ii^)2  =  2500  iiSl!^   , 

(  sec) 

where  Manning's  n^'  .03,  was  suggested  by  the  engj.neers.   The 
differential  equations  are  taken  in  the  form  (see  p.  lii. 
Report  I,  eq.  ( 3.6) ) 


(2.1) 


2cc„   +  V.   +  vv„  +  E  =  0 

X        T>  X 


2c,+2vc   +CV      =0 

t        XX 


where 

E  =  -gS+  gS^  5   V  =  velocity,   y  =  deptn  of  water  , 

c  -  Jf7   (the  propagation  speed  of  small  disturbances)  . 

Tn  carrying  out  the  computation  it  vjas  found  convenient  to  use 
the  mile  as  unit  of  length  and  the  hour  as  unit  of  time. 

Tile  physical  problem  considered  wa.?-  the  following:   at 
time  t  =  C,  steady  flow  of  depth  20  ft.  is  assumed.   At  the 
"head",  X  -  0,  of  the  river,  we  impose  a  linear  increase  of 
depth  which  brings  the  level  to  I].0  ft.  in  1^.  hours.   We  then 
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maintain  the  level  of  bO  ft.  at  x  =  0.   The  initial  velocity  of 
the  water  corresponding  to  a  uniform  flow  of  depth  y  -  20  ft. 
is  calculated  from  ."-,  =  S  to  be 

V   =  2,38  mph   } 

the   propagation   speed   of   small   distiortaances    corresponding  to 
the   depth  of  20   ft.    is 


c 


o 


=  /gy7  =  17.3  mph  . 


The  problem  then  is  to  determine  the  solution  v(x,t),  c(x,t) 
for  all  later  time,  t  >  0,  along  the  river,  x  >  0.   Figures  1 
and  2  present  the  result  of  the  computation  in  the  form  of  stage 
and  discharge  curves  plotted  as  functions  of  distance  along  the 
river  at  various  times.   As  stated  above,  we  vjorked  with 
the  mile  as  the  unit  of  length,  and  the  hour  as  th<;-  unit  of 
time.   Our  stage  measurements  were  then  converted  into  the 
customary  unit,  the  foot, and  the  discharge   into  cubic  feet 
per  second. 

2.2   Compiitational  techniques.   In  order  to  indicate  how 
the  solution  was  calculated  it  is  convenient  to  refer  to  a 
diagram  in  the  (:c,t)  plane  (also  see  FigiAre  3)' 


-■*:.?^.     i 
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1.25 


.7 


hourr 


VI 


miles 


Diagram  1 


■Regions    In  which  various    computational  methods  were   tried 


According  to   the    theory  presented   In  Repoi't   T,   we   know  that 

for  X   >    (v^+c    )t   =  19.7    t,    called  re;- -Ion  0,    the    solution  Is 

given  by   the   unchanged   Initial   data,    v(x,t)    =  v,    c(x5t)    =   c 

(since    the    forerunner   of   the    distiorbance    travels    at   the    speed 

V    H-  c      =19.7  mph)  . 
00 

Tn  order  to  compute  the  solution  by  finite  difference 
methods,  we  found  it  necessary  to  use  rather  small  intervals 
near  the  origin.   This  is  to  be  expected  since  the  values  of 
c.  were  changed  discontlnuously  at  t  =  0,  x  =  0.   We 
experimented  with  various  Interval  sizes  and  finite  difference 
schemes  in  order  to  determine  the  least  time  consuming  way  to 
calculate  the  progress  of  tie  flood.   Our  conclusions  have 
been  summarized  in  the  introduction  to  this  report,  ^e    shall 
now  describe  the  various  schemes  and  the  regions  in  which  they 
were  used  (see  Diagram  1  and  Figure  3). 
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Region   I,      0<x<19.7t,    u<t<    .}\..      Intewals   /\x.  --   1  mile 
and   /\t  -    .  0[|.8  hour  wore    required   owing   to   the    sudden   Incr 


•p  r>  p.  r.  I 


of  depth  at  x  =  0,  t  =  0.   The  finite  difference  formulae 
derived  in  Report  I,  eqs.  (I).. 8),  (1+.9)  were  used: 


Vp  =  Vj^  + 


At 
Ax 


-  (  Si  -  ^M^  ^7  ^M  -  I  ^R  -  Si  ■"  S  ^  ■  ^''^M 


(2.2) 


^P  -  c^  +  •? 


1  At 
Ax 


^M-'VHt  Vj.  -^  v^  +  c^_-c„j) 


■^  (s  -^>i^(^^m4^r-s-^s.) 


where  P  refers  to  a  point  A^  later  than  M,  while  L  and  R  are 
/\x  xinits  to  the  left  and  right  of  M  respectively. 


rt 


■A^ 


X 


At 


Ax- 


R 


x 


Diagram  2 


Net  points  used  in  rectangular  scheme 

(VJe  have  adopted  the  convention  of  using  v  .  to  denote  the  value 
of  the  function  v  at  the  point  M.) 


•/  ••■  -  ,  ■/   ^  K.  ,v 


'".i  :•>  '< 
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In  Region  II,  0<x<19.7t,  ,I|.<t<.7,  with  /\x  =  1  mile, 
A  t=  .  02I4 ,  vje  sxperlmented  with  a  different  finite  difference 
scheme,  called  a  "staggered"  scheme.   This  method  is 
appropriate  for  use  when  the  space  intervals  are  uniform  in 
size  along  the  length  of  the  river.   The  scheme  is  described 
in  Report  I,  eq.  (K.ILl)  and  eliminates  any  reference  to  values 
at  the  point  M. 


At 


L  ♦<—  A  X ^-^ 


Ax  — —  R 


X 


Diagram  3 
Net   points   used   in  "staggered"    scheme 


(2.3) 


where 


At 


2(^L-°R^°M^^^^\,  -^R^^M^^-A^^M^. 


1  At 


-^°L-"R^V^^^^L"^R^°r!^^ 


M^=- 


=   s^^L-^^R^    '      V 


=  |(Cj^-»-Cp^)     ,      E,,,,  =  |(El+  J^^) 


¥e  must  pay  a  slight  penalty  for  using  this  somewhat 
simpler  formula.   That  is,  we  find  that  in  order  to  calculate 
the  value  of  v(0,t),  the  velocity  at  the  head  of  the  river,  we 
must  introduce  some  additional  consideration.   ^\'e  refer  to 


t  p 

■M 


R 


B 


G 
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H 


^ 


Diagram  ii 


Net  point  arranp:ement  used  at  boundary  in  "  ^tan;g:ered"  scheme 


Diagram  kt   where  we  indicate  schematically  the  staggered  net 
points  P,  R,  F,  G  and  H.   It  is  clear  that  equations  (2,3) 
together  with  the  prescribed  values  of  c(0,t)  do  not  produce  a 
formula  for  Vp  in  terms  of  quantities  evaluated  at  the 
preceding  staggered  points  R,  F,  G,  H,  etc.   '7e  can  overcome 
this  difficulty  by  determining  values  for  Vg,  Cg  by  linear 
interpolation"  from  the  values  at  the  points  P  and  G  (or  if 
rapid  changes  are  made  in  the  boundary  valves,  we  may  use 
quadratic' '  interpolation  on  the  values  at  the  three  points 
F,  G,  H) .   That  is,  we  use  such  determined  values  of  v^,  c^  in 


For  linear  interpolation 


For  quadratic  interpolation 


^B  =  2-^°F^^G^   • 


'B 


(3c„4- 6c^  -c„)/8  ,  v_  =  (3v^  +  6v  -V  )/8   , 


G   H 


we  used  these  formulae  to  smooth  the  data  near  the  boundary 
when  we  changed  interval  sizes. 


.')l 
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tha  difference  equation  {2,2)    to  determine  the  value  of  v,^. 
We  then  use  (2,2)  again  to  determine  Vp.   The  necessity  for 
using_a  different  formula  to  calculate  at  boundary  points 
would  add  seme  difficulty  to  a  scheme  for  automatic  computation, 

Region  TIT,  C  <  x  <  5  miles,  .7  <  t  <  1.25  hr,  v;lth 
A  ^    =1  mile,  At=  .02I4.   hr .  We  used  the  same  formulae  as  in 
Region  II  (except  for  points  on  the  boundary  between  Rep-ions 
III  and  I^S  see  below). 

Region  IV,  5  f  x  <  19.7  t,  .7  1  t  <  1.2^,  with  A,x  =  2  miles, 
A  ^  -    .Oli-B  hr .   The  values  at  the  ooundary  between  Regions 
III  and  IV  were  obtained  by  linear  interpolation  from  the 
neighboring  values.   Other  quantities  were  computed  by  the 
"staggered"  scheme  as  in  Regions  IT  and  III. 

Region  V,  0  <  x  <  Vt ,  1.2^  <  t  <  10  hrs,  /\x  =  5  miles, 
At  =  .17  hr .   V  represents  a  variable  speed  which  marks  the 
downstream  end  of  the  observable  disturbance  (V  =  10) .   Ihat  is, 
by  using  an  expansion  scheme  (see  Appendix  I)  we  determined  a 

Region  VI,  Vt  <  X  <  19.7  t,  about  the  forerunner  of  the 
disturbance,  in  which  the  floi-;  is  practically  undisturbed  to 
the  accuracy  we  were  Interested  in.  Vie   therefore  used  this 
expansion  to  calculate  the  various  qi.i.antities  in  Region  ^'I 
and  used  the  staggered  scheme  to  compute  the  values  in 
Region  V.   ^'e  permitted  Regions  V  and  i^I  to  overlap  slightly, 
in  order  that  (a)  '.-je  should  a\''oid  having  to  Interpolate  for 
boundary  values  aad  (b)  that  we  might  have  a  check  on   the 
consistency  of  tha  tvjo  methods  of  computation. 


iJo,:.??^    f'f:: 


?-.,-(%        T 
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2*3   Observations .   The  conclusions  reached  on  the  basis 
of  these  calculations  of  a  flood  in  the  Ohio  can  be  summarized 
as  follows: 

(a)  'Ihe  model  flood,  with  a  5  foot  per  hour  rise,  is  extreme 
even  in  comparison  with  the  191+5  flood  for  which  the  rate 
of  rise  at  i-Fheeling,  West  Virginia  was  less  than  .7  foot 
•oer  aour .   Such  a  case  exaggerates  the  way  in  which  our 
finite  difference  methods  propagate  small  errors.   For 
example,  we  found  slight  inaccuracies  developing  at  the 
head,  x  =  0,  when  we  increased  the  /\x  interval  size. 

Vfe  could  have  devised  special  techniques  to  smooth  this 
effect--but  vie  didn't  in  order  that  we  could  better 
observe  how  the  errors  were  propagated. 

(b)  In  spite  of  the  e'lceptionally  high  rate  of  rise,  the 
inaccuracies  created  by  using  our  finite  difference 
methods  were  damped  out  rather  strongly  (in  about 
8-10  time  steps,).   it  is  possible  to  control  these 
inaccuracies  by  using  appropriate  smoothing  techniques 
or  more  simply  by  using  small  interval  sizes. 

(c)  The  ■orocess  by  which  the  small  errors  of  the  finite 
difference  scheme  die  out  may  be  described  as  follows: 
A  larger  value  of  v  produces  a  larger  friction  force 
V7hich  alowa  down  bhe  motion  and  produces  at  tho  n-oxt  time 
a  smaller  velocity.   The  lower  velocity  in  a  similar  way 
then  operates  through  the  resistance  to  create  a  larger 
velocity  and  the  process  repeats  with  a  steady  decrease 
in  the  amplitude  of  variation. 


■fi  ?!C.':  i^rn*-'-  v- J     i :,:  ;:■ 


^ 
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(d)  'fe  checked  the  accuracy  of  our  computation  (as  a  function 
of  the  interval  size)  oy  repeating  the  calculation  for 
various  Interval  sizes  over  the  same  rep:ion  in  space  and 
time.   Such  self-checking  by  changing  intervals  can  be 
coded  into  an  automatic  procedure  for  computation. 

(e)  The  "staggered"  scheme  introduced  additional  diff iciilties 
for  the  computation  of  boundary  values  (in  contrast  with 
the  rectangular  schemes) .   This  trouble  was  emphasized 

in  the  computation  at  the  boiindary  of  Regions  III  and  IV 
(smaller  A^x  intervals  were  used  in  Region  III).   Our 
plan  to  use  the  stage  gaging  stations  as  net  points  would 
be  simpler  to  code  with  a  rectangular  scheme;  and  would 
not  reaulre  any  special  treatment  at  the  boundaries. 

(f)  Iffe  have  yet  to  investigate  the  possibility  of  using 
"implicit"  finite  difference  methods--which  may  permit 
the  use  of  larger  time  steps--at  the  extra  cost  of  using 
r.ore  complicated  procedures  to  solve  the  equations  at  each 
time  step. 

(g)  The  linearized  theory,  obtained  for  a  small  perturbation 
about  the  uniform  flov/  with  20  ft.  depth,  did  not  give 
an  accurate  dsscrlDtion  of  the  solution  of  our  problem. 
We  compared  the  stage  profiles  rsredicted  by  the  linear 
theory  and  those  calculated  from  the  non-linear  equations. 
After  a  period  of  less  than  ?.   hours,  there  was  a  deviation 
of  2  feet.   The  linear  theory  predicted  smaller  values 
for  stage  and  dj-scharge. 


::C:.'.'  '"i:'-.  r'.ri'i.". .  '■fj'i-  ■■•'-)>/ 
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(h)   It  iNfOuld  be  convenient  to  be  able  to  obtain  in  advance 
a  safe  estimate  for  the  maximum  value  of  the  particle 
velocity,  in  order  to  select  an  appropriate  safe  value 
for  At,  since  we  must  have  At  <  ==-^' (  see  Report  I, 

: — L   '  i — ■■       —  V+C  " 

P.  26).   Our  calculations  ".ndicate  that  thi?  may  not 
be  too  simple  theoretically,  for  the  maximum  velocity 
at  X  =  0  greatly  exceeds  its  asymptotic  value. 

v(0,t) 
I   mph 

i  .■•        X 


/ 


3.7  '--^ 
I  / 

1/ 


i 


velocity  for  \\0   ft, 
steady  flow 


'"  t  hours 


h 


Dia.eram  5 


^'Jater  velocity  obtained  at  "head"  of 
river  thro\;Ah  re^':ulation  of  stage 


But,  it  may  vjell  be  possible,  on  the  basis  of  empirical 
observations,  to  obtain  adeouate  estimates  of  the  maxim\Am 
velocity.   As  was  indicated  earlier,  xje  can  easily  code 
the  calculation  of  the  proper  interval  size  as  a  part 
of  the  computation  should  other  means  of  deterxi'ination 
be  inadequate. 
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(i)   ■'■'e  note  that  the  curves  of  constant  stage  have  slopes 

which  are  closer  to  5  J^-Ph  (the  speed  iilth  which  a  steady 
progressing  flow,  1;0  ft.  upstream  and  20  ft.  downstream, 
moves,  see  next  section)  than  to  the  over  19  mph  speed 
of  propagation  of  small  disturbances. 


I 


t  hours  ..   re/xion   of 

/    practically 


/ 


10  ' 


undisturbed 

/  ■'■     /    flow  of  20' 

/    /    ■'■     .  ■     / 
I 

/    /     ..'     / 
.'pp.'  ./  ^M-  '  •'  cin'      --""^  ■ 


5  mph- slope     19.7  mph- si  ope 

19.7  -^\--[Wl 


^0  100  X  miles 

Diagram  6 

Curves  of  constant  stage--comparison  X'jxth  f i rs_t 
characteristic  and  steady  progressing  flow  velocity 


The  region  of  practically  undisturbed  flow  (determined 
by  expansion  about  the  "first"  characteristic  x  =  19.7  t, 
seo  Appendix  I)  is  shown  above.   In  practice,  v/e  expect 
the  local  runoff  discharges  and  the  non-unifcrm  flow 
conditions  to  eliminate  the  region  of  practicalljr 
undisturbed  flow.   For  this  reason,  in  general,  we  can 
not  use  analytic  expansion  schr-ijes  to  tr;'"  to  save 
computational  labor  of  this  kind. 
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2,[\.     Determination  of  accuracy  of  finite  difference  method 
by  comparison  with  an  exact  solution.   Ihere  is  a  case  in 
which  an  exact  solution  of  the  differential  equations  is 
knovjn,  i.e.  the  case  of  a  steady  progressing  wave  with  two 
different  depths  at  great  distances  upstream  and  downstream. 
This  exact  solution  for  the  case  of  a  v/ave  of  depth  20  ft. 
far  dot-mstream  and  liO  ft.  fi^r  upstream,  was  taken  as  furnishing 
the  initial  conditions  at  t  =  Of  or  a  vjave  motion  in  the 
river.   Mith  the  initial  conditions  prescribed  in  this  vray 
the  finite  difference  method  was  used  to  determine  the  motion 
at  later  times;  of  course  the  calculation,  if  accurate,  should 
furnish  a  wave  profile  and  velocity  distribution  which  is  the 
same  at  time  t  as  at  the  initial  instant  t  =  0  except  that 
all  quantities  are  displaced  downstream  a  distance  IJt,  with 
U  the  speed  of  the  steady  progressing  wave.   In  this  way  an 
opportunity  arises  to  compare  the  approximate  solution  with 
an  exact  solution.   In  the  present  case  the  phase  velocity  U 
is  ar) proximately  5  mph.  ^-^e  considered  an  interval  size  of 
A.X  =  10  miles  in  a  "staggered'*  scheme  viith  /\t  -    .31^  hr . 
After  7  hours,  the  calculated  stage  values  a^'-eed  to  vjithln 
.6  per  cent  with  the  exact  values.   The  dischar^^e  and  the 
velocity  deviated  by  less  than  2.5  per  cent  from  the  exact 
value.   The  fact  that  stage  values  can  be  predicted  more 


In  Appendix  II  we  derive  the  formulae  which  describe  the 
steady  progressing  wave. 
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accurately  than  discharge  values  is  a  characteristic  of  the 
presently  used  finite  difference  iiethods  that  was  observed 

in  all  calculations. 

!■  y  feet 

ko 

20  - 


I 


X 


3.7t 


,  V  mph  100  miles 


2J4r 


'  Diagram  7 

Steady  progressing  flow--stage  and  velocity  profiles 


^*   Ihe  junction  of  the  Ohio  and  the  Mississippi 

3.1   formulation  of  the  junction  problem.  Oiir  model  for  the 
junction  of  the  Ohio  and  Mississippi  Hivers  was  constructed 
to  test  our  finite  difference  methods  with  regard  to  their 
applicability  in  such  cases.   We  believe  that  the  consideration 
of  the  actual  problem,  as  it  occurs  in  nature,  would  present 
no  new  difficulties  in  principle  (except  for  the  consideration 
of  overbank  flow).  We  supposed  the  upstream  Mississippi 
section  to  be  identical  with  the  Ohio  River — rectangular  cross- 
section,  1000  ft.  wide,  slope  of  .5  ft. /mile,  Manning's  n=.03, 
infinite  in  length.   We  imagined  the  downstream  Mississippi 
branch  to  be  rectangular,  twice  as  wide,  i.e.  20"0  ft., 
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Manninp:'s  n=.03,  and  a  slightly  smaller  slope  .k9   ft. /mile. 
The  latter  modification  of  the  slope  was  made  in  order  to 
obtain  an  initial  solution  corresponding  to  a  uniform  flovj  of 
20  ft.  depth  in  all  three  branches  (such  a  change  is  necessary 
to  overcome  the  decrease  in  vretted  perimeter  which  occurs  on 
going  downstream  through  the  junction) . 

[  2  ]    \ 

TTp  stream  \ 
\  Mississippi 

\  \ 


-  1000' 


Downstream 
Mississippi 


Diagram  3 
Schematic  plan  of  junction 


I'fe  considered  the  problem  of  starting  a  flood  wave  in  the 
Ohio,  in  fact,  the  same  flood  we  had  computed  earlier.   The 
difference  was  that  we  began  the  Ohio  flood  $0   miles  uostream 
from  the  junction.   After  about  2.5  hours  the  forerunners  of 
the  Ohio  disturbance  reached  the  junction.   They  then 
separated  into  two  wavelets,  one  going  upstream  and  the  other 
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downstream  on  the  Mississippi.   In  addition,  at  the  same 
instant,  a  reflected  wavelet  started  backwards  up  the  Ohio. 
Our  finite  difference  calculations  were  begun  in  all  three 
branches  from  the  moment  that  the  junction  was  I'eached  by 
the  forerunner  of  the  Ohio  flood.   ('Je  saved  the  labor  of 
recomputing  the  motion  for  the  first  2.5  hours  in  the  Ohio, 
since  the  motion  was  identical  with  the  one  previously 
calculated.)   The  flow  was  followed  for  10  hours  and  the 
results  are  indicated  in  the  stage  and  discharge  plots  as 
a  function  of  distance  (for  various  times)  given  in  Figures 
[}  and  5  • 

3 . 2  Description  of  the  computational  procedur e s .   Let 

^ni'  ^11]'    ^\2]'    ^12]'    "^[31'  ^[■^l  ^spresent  the  velocity  v 
and  propagation  speed  c  for  the  Ohio,  upstream  Mississippi, 
and  downstream  Mississippi,  respectively.   A  "staggered" 
scheme  was  used  with  interval  /\^x    =5  miles  and  /\^t     =  .17  I'lr 
as  indicated  in  Diagram  9.   The  junction  point  is  denoted  by 
X  =  0,  the  region  of  the  Ohio  and  the  upstream  Mississippi 
are  represented  by  x  <  0,  while  the  downstream  Mississippi 
is  given  on  x  >  0.   The  time  t  =  2.5  hr .  corresponds  to  the 
instant  that  the  forerunner  of  the  flood  reaches  the  junction. 
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Diagram  9 
Junction  net  point  scheme 


The  values  of  the  Quantities  v  and  c  at  the  junction 
were  determined  as  follovrs:   Let  us  assume  that  the  values 
of  V  and  c  have  been  obtained  at  all  net  points  for  times 
preceding  that  of  the  boundary  net  point  P.   We  knovj  that  at 
the  junction 

"=[1]  "  "^[2]  "  "^[3] 

(since  c  =  Jgy     and  the  three  branches  are  level  at  the 
junction)  and 

y[l]^[l]  ^  y[2]^[2]  =  ^^[3f[3l 


■•;>i*. 
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(since  what  flows  into  the  .iunction  from  the  upstream 
Mississippi  and  the  Ohio  must  flow  out  of  the  junction  into 
the  downstream  Mississippi).   If  the  values  of  v  and  c  were 
known  at  M  in  the  respective  branches  of  the  rivers,  we  could 
find  the  values  at  P  from  the  following  three  equations: 
Report  I,  eq.  (1|.6)  for  the  Ohio  and  upstream  Mississippi, 
and  Report  I,  eq.  (1+.7)  for  the  doimstream  Mississippi.  We 
rewrite  all  of  the  equations  for  convenience 


Cp[i]  =  Cp^2]  =  °F[3]      (with  c  =  v/gy  ) 

^P[l]    ■"   ^^P[2]    =  2Vp.,^  (sine,   y^i^    =  j:^^    =  j^^^) 


(3.1^2]^        ^ 


I  ^F[2]-°M[21    ,    .  ^„         )    ^'^MLajLAlljl) 

I AT ^'n[2]'M2]^         '^,  V 

j'^-I2llZli[21    ,    ,,  ^^  )  r_ML2lllLL2I^\     .  E  =0 

I  A^  ^    ^^U2]-^M[2])^  Ax  J]  ^f25  ' 


and 
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The  above  system  of  six  equationa  determines  uniquely  the 

value  s  V  r  -j_  -J ,  c  ^  -j_  -. ,  v  j-  2  j ,  ^  [  p  ] '  ^'  [  3  ] '  ^-  [  3  ]  ^'^'^  ^  -^^  terns  of 
their  values  at  the  preceding  point?  L,  ''I  and  p.   -The  equations 
can  be  solved  explicitly.   On  the  other  hand,  th'^  values  at  M 
are  determined  in  the  same  way  from  the  preceding  values  at 
A,  F  and  B.   The  values  at  A  and  B  are  determined  by  linear" 
interpolation  between  the  nel.^boring  points  (K,P)  and  {P,G) 
respectively  (see  Diagram  9). 


3 . 3   Observations . 
(a)   I'e  find  that  the  flow  tend;?,  to  a  steady  state  in  which 

the  downstream  •'•'Mississippi  has  uniform  flovj,  as  indicated 
below  (see  Appendix  TI). 
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Asymptotic  flow  in  junction  problem 


Ouadratic  interpolation  using  (F,G,H)  and  (J,K,F)  was  used 
after  the  stage  increased  appreciably  at  the  j\mction. 
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(b)  The  existence  of  both  an  up- and  downstream  branch  of  the 
river  into  which  the  Ohio  discharges  has  the  effect  of 
lowerlnp  the  Ohio  River  level.   A  crude  estimate  of  what 
would  happen  should  the  Ohio  discharge  only  into  the 
dov;nstream  branch  of  the  Mississippi,  indicates  a  rise 
of  about  5  per  cent  in  the  level  of  the  Ohio  at  the 
junction. 

(c)  There  is  a  definite  influence  on  the  upstream  branch  of 
the  Mississippi  (see  Figures  [j_  and  5  )  and  it  is  in  this 
connection  that  the  finite  difference  methods  display 
their  advantage  over  the  standard  flood  routing 
procedures.   For  example,  at  a  point  $0   miles  upstream  in 
the  Mississippi  the  stage  is  increased  by  1  ft.  as  a 
result  of  the  pouring  of  the  Ohio  River  flood  into  the 
junction. 

(d)  i''e  found  it  convenient  to  compute  the  flow  in  the 
neighborhood  of  the  flood  forerunners  in  the  two  branches 
of  the  Mississippi  by  means  of  the  expansion  scheme 
mentioned  earlier  and  described  in  detail  in  Appendix  T. 

(e)  The  damping  of  small  errors  and  other  features  peculiar 
to  these  equations  are  the  same  in  the  present  case  as 
in  the  preceding  case  of  the  Ohio;  they  are  described  in 
section  2,3  (in  connection  with  the  Ohio  River  model 
flood  calculation). 


.I.*.  I 


>;V'V.     'V! 


j.:-f-t     <■ 


On.:"'    ^i'    olr'^j    '^o-'    '"io     C^V"^ 


^;.'r:l  .7i 


■^    NrH^    '  ^  ^':-:^    r 


'■■■         "    '.   !  ^     _         ■- 

I-    vi.iJ.C'Cr    to   'J  I-    , 


J . 


:vo   0  7 


'  • '    e 


33 


o 
c 


a> 


CJ» 


o 

\_ 
o 


o 

Q. 
O 


T3 

C 

a 
o 

O 


en 


-    c 


a» 

E 

o 

X- 

CO 

a» 

.»_ 

if) 

3 

<u 

c 

o 

n 

O 

-J 

"  — 

c 

a* 

J= 

<v 

3 

E 

XJ 

u 

—> 

c 

o 

c 
o 

to 

E 
o 

0; 

CP 

O 

F 

o 

<4— 

\_ 

^— 

»♦— 

TD 

M— 

to 

■♦— 

o 

II 

M 

II 

-T    ~       -J^— -<-'«-  -"^--4—  —-    •      "Ol         •- 


t- 


o 


Q. 
Q. 

(/> 
</) 


4- 


r 


i-^ 


'.,T 


o 


o 


O 
CVJ 


O 


c 
o 


34- 


c 


35. 


Appendix  I 

Expansion  In  the  neighborhood  of  the  first  characteristic . 

It  has  been  mentioned  already  that,  whereas  the  first 
signal  of  a  disturbance  initiated  at  a  certain  point  in  a 
river  with  uniform  flow  travels  downstream  with  the  speed 
V  +  Jgy^,  the  main  part  of  the  flood  wave  travels  more  slowly 
(cf.  D^ymie'),  depending  strongly  on  the  resistance  of  tne 
river  bed.   An  investigation  of  the  motion  near  the  head  of 
the  wave,  i.e.  near  the  first  characteristic  with  the  equation 
X  =  ^(v  +c  ),  shox'js  immediately  why  the  main  part  of  the 
disturbance  will  in  general,  but  not  always,  fall  back  behind 
the  forerunner • 

TLae  motion  is  investigated  in  this  Appendix  by  means 
of  an  expansion  of  G.  vrhitham  anc"  A.  Troesch  and  is  carried 
out  to  the  two  first  orders  for  the  model  of  the  Ohio  River, 
and  to  the  lowest  order  in  the  much  more  complicated  case 
of  the  junction  problem.   The  results  obtained  enable  us  to 
improve  the  accurac?/  of  the  solution  near  the  first 
characteristic.   It  turns  out  that  the  finite  difference 
scheme  yields  river  depths  which  are  too  big  (see  below), 
owing  to  the  fact  that  om^  large  mesh  width  smoothed  the 
river  profile  out. 


Deymie,  Fh.  Propagation  d'une  intumescence  allonges  (problems 
aval).   Int.  Cong,  for  Appl.  Fech.,  5th,  Proceedings, 
Sept.  1938,  N.  Y".  ,  ^^'iley,  1939. 
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Diagram  11 
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Error  introduced  by  finite  difference  scheme 
in  neip:hbor'hood  of  first  characteristic  of  a 
rapidly  rising  flood  wave 


In  order  to  expand  the  solution  in  the  neighborhood  of 
the  wave  front,  we  introduce  new  coordinates 


=  x   and   t:  =  (v  t- c  )t  -x 


such  that  the  r-axis  (i.e.  Z~   0)  coincides  with  the  first 
characteristic.   The  basic  system  of  equations  (?.l)  are 
restated  for  convenience 


(A.l) 


2cc  +v.  +  vv^,  -  tS  +<r,Sf.   =  0   , 


cv^  +  2vc^^  +  2c^ 


=  0   . 


Upon  substitution  of  the  new  variables  ^  and  t    we  find 
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2c(c^-  c^)  +v(v^  -  v^)  +     (Vq+  c^)v^- gS +gS^  =   0      , 


where    the    friction    slope  S^for   a  rectangular   channel   of  width 
B   is   given  by 

vivl  g^^/3      ,    1/12  i^^/3 

S  ,     = ; — —    =     — V    V  M    — g-  +  — —  > 

•      X-V^  ^  1  c^     S3  ( 

1  + 


B 


'.'e  expand  v  and  c  as  power  series  in  x    with  coefficients  that 
are  functions  of  5  as  follows: 

V  =  v^  +  v^ ( ^ ) r  +  V2( - )r^  +  . . .   , 

2 

c  =  c   +  c^(  t )  1:;  +  C2(  c  )'i.'  +  .  .  . 

This  expansion  is  to  be  used  for  f- -•  0  only,  since  Tor  tf  <  0 
we  are  in  the  undisturbed  region  and  all  fu.nctions 
V]^(.^,)  ,V2(C)  >  •  •  •  J  c^(c;j  ,c^(^)  ,  .  .  .   vanish  identically.   If  we 
insert  the  series  for  v  and  c  into  equations  (A.l)  and  collect 
terms  of  the  same  order  in-  T,  we  get  ordinary  differential 
equations  for  v-,(fc),  c,(f),...   for  example 

^°^      -         ^  -  2 L—l  = 


. .   ,  --  0 


MVo-o)'-dl^-^2c^Uc^.gslJ--3^^    ^^^  ^ 

^  ""'  gB 


Although  the  solution  of  this  differential  equation  for  c^(£;) 
is  easily  obtained,  the  result  expressed  in  general  terms 
looks  so  complicated  that  we  prefer  to  give  it  only  for  the 
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case  of  the  Ohio  River  under  the  conditions  discussed  in 
section  2:   c,  =  ( 1 ,  05  +  3 .  06  e"*'  ^  ^)~  $    in  niles  and  hours. 
This  result  has  the  following  physical  meaning:   The  angle  a 
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Diagram  12 


of  the  profiles  measured  between  the  vrav-?  front  and  the 
undisturbed  water  surface  is  shown  to  div  out  exponentially 

a>w ^  ,  with  a  and  b  constants  dep-endinj;^  en  the  river 

1  +  ae 
and  the  boundary  condition  at  x  =  0,   Theoretically,  a  could 

also  increase  exponentially  downstream  so  thr?t  a  bore  v.'ould 
eventually  develop,  but  only  if  the  increcr-o  in  level  at  x  =  C 
is  extremely  fast;  in  our  example  (sec.  2,1)  no  bore  will 
develop  unless  the  ^^'ater  rises  at  a  rate  of  at  least  1  ft, 
per  minute. 

Unfortunately,  the  evalucction  of  C2(C)  which  yields 
the  curvature  of  the  profile  at  the  v;ave  front  is  already 
very  ciimbersome.   The  curvature  is  found  to  decrease  for 
large  x  like  xe'   ,  b  being  a  positive  constant.   With  the 
two  hip;hest  order  terms  in  the  expansion  Jcnown,  we  were  able 
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Diagram  13 

Reg;ion  of  practically 
undisturbed  flow 
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mile! 


to  estimate  the  region 
adjacent  to  the  first 
characteristic  where  the 
flow  is  practically- 
undisturbed.   It  is 
remarkable  how  far  behind 
the  theoretical  forerunner 
the  first  measurable 
signal  travels,  see 
Diagram  13  • 

Tn  a  similar  way,  an 
expansion  as  a  power  series 
in  K     has  been  carried  out 


for  the  junction  of  the  Ohio  and  Mississippi,  as  described  in 
section  3.   Here  even  the  lowest  order  term  required  a 
complicated  computation,  since  vje  had  to  work  simultaneously 
in  three  different  x,t-planes,  with  boundary  conditions  at 
the  junction.   'Ihe  differential  equations  for  c^  are,  in  all 
three  branches,  of  the  same  type  as  for  the  Ohio,  and  their 
solution  for  the  junction  problem  as  treated  in  section  3  ai^s 
C-,  =  .OOOSIj.  e^^^^  for  the  Mississippi  upstream, 
c-  =  .OOO6I4.  e"'"  "^^  for  the  Mississippi  downstream,  in  miles 
and  hours.  This  means  that  the  angle  a  dies  out  exponentially 
again  in  the  "  Mississippi,  a  little  faster  downstream  than 
upstream,  as  could  be  expected,  since  the  oncoming  water  in 
the  upstream  branch  makes  the  wave  front  steeper. 
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In  the  problem  of  the  idealized  Ohio  Hiver  and  of  the 
junction  the  expansions  wer-e  carried  out  numerically  in 
full  detail  and  were  used  to  avoid  computation  by  finite 
differences  in  a  region  of  pi-'actically  imdisturbed  flow." 
At  the  same  time  the  accuracy  in  the  region  of  small 
disturbance  was  improved  in  general  (as  indicated  above). 


This  would  become  more  and  more   important   if  the  flew 
were  to  be  computed  beyond  10  hours.  '  . 
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Appendix  II 

Steady  floijs  and _  steady  progressing;  waves . 

II, 1   Steady  flows .   'fe  define  a  steady  flow  to  be  one  ^or 

which  the  velocity  v  and  depth  y  (or,  e qui va lent Ij'-  the 

propagation  speed  c  =  v/'gy  )  are  independent  of  the  time,  that 

is,   V,  =  y,  =  0.   It  follows  from  the  equation  of  continuity 
(see  eq.  (3»1)  in  Report  I), 

^t  ""  ^^x  +  y^x  =  °   ' 

that  for  steady  flow 

(II.l)    (vy)   =  0  whence   vy  =  D   (Da  constant)   . 
Similarly,  the  eouation  of  motion  (eq.  (3»2)  of  Report  I) 
V,  +  vv^  +  gy  +g(>'?f--S)  =0    yields 

(TI.2)       ^'Vx'^S^x  +g(Sf  -S)  =  0 


Since  from  equation  (II.l)  it  follows  that 


V  =  I   and   v^  =  -  4  y^   , 


equation  (11,2)  becomes 

.2 
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(II. 3)      g-^ly^  -Hgl -^ j^  -S  I  =  0 
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v^e  note  that  eouation  (II. 3)  has  the  simple  solution 
y  =  constant,  for  y  satisfying 


r:  ^ 


,7  •  {.   / 


■        V  :;0-' 


l;2. 


T.2    o  /   ,,   N  U/3 
\  X  t-  g  / 

Ihls  means  that  we  can  find  a  flow  of  uniform  depth  and 
velocity  having  a  given  discharge  BD  (B  is  width  of  channel). 
Conversely,  by  fixing  the  depth  y  we  can  find  the  discharge 
from  (II. I)-)  of  the  corresponding  uniform  flow. 

For  a  channel  with  varying  phj'-sical  parameters  such  as 
cross-section,  friction  term,  etc.  the  steady  flows  provide 
the  ^^Jell-known  backwater  curves.   In  general,  v;e  would  find 
non-trivial  steady  solutions  y  =  y(x)  and  v  =  v(x)  for  the 
non-uniform  channel.   The  explicit  determination  of  the  stage 
and  discharge  would  be  possible  by  numerical  integration  of 
ordinary  differential  equations. 

The  analysis  of  the  steady  flox-j  through  the  junction  of 
the  Ohio  and  Mississippi  Rivers  was  made  for  the  conditions 
cescribed  in  section  3.   That  is,  at  $0   miles  upstream  in  the 
Ohio  the  stage  was  held  at  liO   ft,  while  far  upstream  in  the 
Mississippi  the  stage  was  kept  at  20  feet.   The  latter 
condition  implies  that  the  slope  y   should  vanish  upstream 
in  the  Mississippi,  or  that  y  and  D  =  vy  satisfy  equation 
(II. 1-1-).   As  before,  let  us  use  subscripts  [1],  [2],  and  [3] 
to  represent  quantities  evaluated  in  the  Ohio,  upstream 
Mississippi  and  downstream  Mississippi  respectively.   'Jhe 
conditions  that  must  be  satisfied  at  the  junction  can  be 
simply  written 

(11. 5)  Yfij  =  Yrg]  =  ^[3]    (=  ^j  ^^  u-i^nown) 

(11. 6)  Dj^^  +  Dj2]  =  2D^^^  ,   for  X  =  0   . 
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Equation  (II. 3)  can  now  be  integrated  in  each  of  the  three 
river  branches  (from  the  junction  point  as  origin). 
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¥e    observe  that  y  as  well  as  Dr-,-,  and  Dr,.-,  must  be  determined 
in  order  that  (II. 7)  can  be  used.   (Drpi  as  observed  above  is 
known  from  the  upstream  condition  in  the  Mississippi.) 

The  flow  in  the  downstream  MisGlssippi  is  easily  seen 
to  be  uniform  (i.e.  depth  Y  =  Jj   and  velocity  are  constant) 
and  this  means  that  equation  (II. I)-)  (with  the  use  of  (II.6)) 
holds  in  the  form  ( if  3  =  1000  ft.) 


¥e  schematically  plot  the 

integrand,  I(y),  of  equation 

(II. 7)  (for  given  D  and 

other  constants)  in  the 

region  of  sub-critical  flo^ir, 

that  is,  where  gy  is  greater 

2 
than  V  .   TJe  see  from 

Diagram  II4.  that 


X  = 
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I(y)dy 


i(y) 


p  /  3 

g  -  D"/V_ 


1  \ 
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will  become  positive  infinite  Diagram  iL). 
only  If  y  =  y  =  value  for 

which  E  =  0.   (To  other  combination  of  finite  values  of  y 

and  y  can  produce  a  positive  infinite  value  for  x.) 
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Equation  (II. 8)  together  with  equation  (II. 7)  rewritten 
for  the  ^Q   mile  segment  of  the  Ohio  as 


yr.-.      ._^il 
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are  a  pair  of  simultaneous  equations  from  waich  the  values 

of  y  and  Dr,-,  can  be  determined  by  an  iterative  procedure, 
J       L  -1-  J 

The  results  are  for  x  =  0, 

^[1]  =  ^[2]  =  ^[3]  =  ^J  ^  ^^'^    "• 

Vr,-,  =  I4-.83  mph  ,   Vrpi  =2.33  mph  ,   Vr^-,  =  3«l8  '^ph 


The  curves  of  stage  versus  distance  were  computed  from 
formula  (II. 7)  and  are  drawn  below. 
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Diagram  15 

Steady  flow  profile  in  neighborhood_of 
junction  of  the  Ohio  and  Mississippi  Rivers 
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I'^e  have  noted  earlier  that  the  transient  solution  of 
our  junction  problem  tends  to  the  above  drawn  steady  flow 
as  time  increases. 

II. 2   Steady-progressing  waves.   For  the  case  of  the  uniform 
channel  it  is  possible  to  find  solutions  in  the  form  of  steady- 
progressing  waves,  that  is 

y(x,t)  =  y(x  -IJt)  ,   v(x,t)  =v(x-Tjt) 

with  U  a  constant  speed  of  propagation. 

In  order  to  find  a  convenient  representation  of  such 
phenomena,  we  introduce  a  new  coordinate  system  (s,t)  by 
setting 

z  =  X  -  Ut  ,     t  =  t 

i^'e  note  that  z  =  constant  corresponds  to  a  motion  with  velocity 
U  in  the  (x,t)  plane.   The  differential  equations  become 

(11.10)  y^  +  (V  -Ujy^  +  yv^  =  0 

(11.11)  v^  +  (v-U)v^  +  gy^  +  g(S^-S)  =  0 

The  requirement  that  y  and  v  be  functions  of  z  only, 
means  that  y,  =  v   =  0  and  as  in  the  case  of  steady  flow  we 
find  that 

(v-U)y=D    (Da  constant)   and 

(TI.12) 
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Eauation  (11.12)  can  be  integrated  to  yield 
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The  integrand  I(y)  for  sub-critical  flows  will  be  typified 
in  the  follo^^Jing  sketch. 
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Diagram  l6 


We  observe  that  there  are  two  vertical  asymr^totes  y  =  yQ>y-| 
between  v?hich  I(y)  is  negative.   This  means  that  for  the 
infinitely  long  river  we  can  obtain  a  non-uniform  steady- 
progressing  wave  by  letting  y.  correspond  to  the  depth  at 
X  =  -co  and  y   to  the  depth  at  x  :=  +co  .   If  D  and  U  are 
prescribed,  y  and  y,  are  uniquely  determined  and  vice  versa 
if  y  and  y,  are  prescribed  then  tt  and  I")  are  uniquely 
determined.   The  latter  viewpoint  is  easier  to  adopt,  that 
is,  we  find  by  setting  E  =  0  for  y  =  y„  and  y  =  y^: 
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Equations  (II.lL|.)  can  easily  be  made  linear  to  determine  IT 
and  D  simply.   In  the  case  we  treated  numericallj^  y-,  =  Lj.O  ft, 
y  =  20  ft,  we  found  that  U  =  5  mph. 

Unfortunately,  the  steady-progressing  waves  do  not  exist 
in  rivers,  since  the  channels  aren't  uniform.   Nevertheless, 
the  soeed  of  propagation,  U,  is  very  close  to  the  observed 
speed  of  travel  of  flood  waves.   For  this  reason,  some 
theoretical  investigations  are  being  made  to  establish  the 
relationship  between  propagation  speeds  of  transient  phenomena 
and  IT,   These  resTilts  are  needed  to  specifj^  the  accuracy  of 
long  range  forecasts  of  flood  conditions. 
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